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Abstract 

The spin- 1/2 Ising model with a spin-phonon coupling on decorated planar lattices partially amenable to lattice vi- 
brations is examined within the framework of the generalized decoration-iteration transformation and the harmonic 
approximation. It is shown that the magnetoelastic coupling gives rise to an effective antiferromagnetic next-nearest- 
neighbour interaction, which competes with the nearest-neighbour interaction and is responsible for a frustration of 
the decorating spins. The strong enough spin-phonon coupling consequently leads to an appearance of the striking 
partially ordered and partially disordered phase, where a perfect antiferromagnetic alignment of the nodal spins is 
accompanied with a complete disorder of the decorating spins. The diversity in temperature dependences of the total 
specific heat is investigated in connection with the particular behaviour of its magnetic and lattice contribution. 
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1. Introduction 

Critical properties of compressible spin systems are tra- 
ditionally subject of extensive theoretical studies, since 
they depend in a rather crucial and subtle way on spe- 
cific constraints laid on a given compressible spin model. 
A vast number of theoretical studies on the compress- 
ible spin- 1/2 Ising planar models serves in evidence that 
the magnetoelastic coupling may either change a contin- 
uous (second-order) phase transition to a discontinuous 
(first-order) one lUl S S B S S lid, cause a renormal- 
ization of critical exponents in a continuous phase transi- 
tion 11211 . or lead to a striking cross-over be- 
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tween both aforedescribed scenarios at a tricritical point 
|ll[ll[ll[ll[ll[ll[ll|23]. From this point of view, 
it is quite valuable to investigate this rather subtle inter- 
relationship between the critical behaviour and specific 
contraints laid on compressible spin systems especially 
with the help of exactly tractable spin models. 

In the present Letter, we will introduce a new class of 
exactly soluble Ising models on decorated planar lattices 
(see Ref. |21] and refernces cited therein) in which atoms 
from nodal lattice sites are excluded from lattice vibra- 
tions in contrast with decorating atoms that are relaxed 
from a condition of perfect lattice rigidity. The lattice vi- 
brations of the decorating atoms will be treated as quan- 
tum harmonic oscillations and will consequently lead to 
a presence of the competing next-nearest-neighbour in- 
teraction, which basically affects the critical behaviour of 
the investigated spin system whenever a relative strength 
of the magnetoelastic interaction exceeds some threshold 
value. However, our assumption on a lattice rigidity of the 
nodal atoms is suflicient to prevent a change in the char- 
acter of a phase transition, which always turns out to be 
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from the standard Ising universality class. 

This Letter is so organized. The investigated model is 
defined in Section |2l where the basic steps of the calcu- 
lation procedure are also explained. The most interesting 
results for the ground-state and finite-temperature phase 
diagrams, the spontaneous magnetization and the specific 
heat are presented in Section[3] Finally, some conclusions 
and future outlooks are presented in Section]!] 

2. Ising model and its exact solution 

Let us begin by considering the spin- 1/2 Ising model 
on decorated planar lattices such as a decorated square 
lattice schematically depicted in Fig. JT] The decorated 
lattice constitute two inequivalent lattice sites to be fur- 
ther referred to as nodal and decorating sites, respectively, 
which are occupied by the nodal atoms with the Ising spin 
cr - 1 12 and the decorating atoms with the Ising spin 
^ - 1 12. Suppose furthermore that the nodal atoms are 
placed at rigid lattice positions contrary to the decorating 
atoms, which are relaxed from the condition of a perfect 
rigidity. Under these circumstances, it is advisable to de- 
fine the total Hamiltonian of the investigated model sys- 
tem as a sum of bond Hamiltonians, whereas each bond 
Hamiltonian "Kt contains all the interaction terms of the 
klh decorating atom and is further divided into the mag- 
netoelastic part 'Hk,m and the pure elastic part 'Tl^jt ; 

Nql2 Nql2 

= £ "H* = J] mk,,n + "Hkj). (1) 

k=\ k=\ 

Above, denotes the total number of the nodal atoms and 
q determines the number of their nearest neighbours. The 
magnetoelastic part of the bond Hamiltonian i^k.m takes 
into account the exchange interaction, which depends on 
an instantaneous distance between the decorating spin /i^ 
and its two nearest-neighbour nodal spins cri-i and 0-^2 

"Hk,,,, = -{J - Aii:)//^(t|j - (7 + Aii:)//^(t|2- (2) 

Here, jl\ and cr^ denote the z-component of the standard 
spin- 1/2 operator with two respective eigenvalues ±1/2 
and the coordinate operator is assigned to a local dis- 
placement of the kt\\ decorating atom from its equilibrium 
lattice position placed in the middle in between its two 




Figure 1: The part of a decorated square lattice. Blue circles label rigid 
lattice positions of the nodal atoms with the Ising spin cr = 1/2, while the 
red ones denote lattice positions of the decorating atoms with the Ising 
spin n = 1/2 that are relaxed from the condition of a perfect rigidity. 

nearest-neighbour nodal atoms. Note that the displace- 
ment Xk will be hereafter regarded to be dimensionless 
(i.e. this parameter determines displacement in multipli- 
ers of some characteristic length), which consequently al- 
lows one to consider all the other Hamiltonian parameters 
to be in energy units. The exchange constant J determines 
a magnitude of the nearest-neighbour interaction between 
decorating and nodal spins provided that the decorating 
atom takes its equilibrium position, while the magnetoe- 
lastic (spin-phonon) coupling constant A describes an ap- 
propriate decrease (increase) in the nearest-neighbour ex- 
change interaction invoked by the elongation (contrac- 
tion) of a relevant distance. The elastic part of the bond 
Hamiltonian 'Hkj includes the kinetic energy of the ^th 
decorating atom with the mass M and the elastic energy 
penalty, which is associated with the displacement of the 
^th decorating atom from its equilibrium position 



It is noteworthy that the elastic part of bond Hamiltonian 
given by Eq. Q is in fact the Hamiltonian of the quantum 
harmonic oscillator with twice as large bare elastic (spring 
stiffness) constant due to a simultaneous distortion of two 
harmonic springs attached to each decorating atom. 

It should be noted here that an interconnection between 
spin and lattice degrees of freedom is usually the main ob- 
stacle, which precludes an exact treatment of the lattice- 
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statistical spin models including the spin-phonon interac- 
tion. However, our assumption on the incapability of the 
nodal atoms to perform lattice vibrations enables one to 
decouple spin and lattice degrees of freedom by making 
use of the local canonical coordinate transformation [2Z] 



Xk 



(4) 



As a matter of fact, the coordinate transformation (|4| 
eliminates from the magnetoelastic part of the bond 
Hamiltonian (|2]l dependence on a displacement operator 



which is substituted through the effective next-nearest- 
neighbour interaction between the nodal Ising spins and 
the less important constant term. A presence of the 
antiferromagnetic next-nearest-neighbour interaction be- 
tween the nodal spins in the effective spin Hamiltonian 
(|5]i implies an eventual spin frustration, which originates 
from the magnetoelastic coupling. By contrast, the elas- 
tic part of the bond Hamiltonian Q remains unchanged 
under the local canonical coordinate transformation dU 



2M 



(6) 



Introducing standard annihilation and creation bosonic 
operators obeying the commutation rule [bk, bl] - 1 



allows one to bring the elastic part of the bond Hamilto- 
nian (|6]l into the diagonal form 



9{:,^hco\blbk + 



1 



(7) 



which is more appropriate for further calculations. It is 
worthy of notice that the angular frequency of normal- 
mode oscillations is given hy oj - y/lK/M. 

Now, let us proceed to a calculation of the partition 
function. It can be easily understood that the total parti- 
tion function can be partially factorized into a product of 
bond partition functions, because different bond Hamil- 
tonians obviously commute with each other. In addition. 



the coordinate transformation (|4|l allows one to treat the 
elastic and magnetic part of the bond partition function 
independently of each other and consequently, the parti- 
tion function can be factorized into the following product 



Nq/2 



Tr,,,exp(-/?K;,,) J] exp(-^; „,) 

K=±l/2 



Nq/2 

^ Y\ -^kjXk,n 
(o-l k=\ 



(8) 



Here, (3 - l/ik^T), is Boltzmann's constant, T is the 
absolute temperature, the symbol ^j^j marks a summa- 
tion over all possible configurations of the nodal Ising 
spins and Tr^ ; denotes a trace over the lattice degrees of 
freedom connected to the kth decorating atom. The elas- 
tic part of the bond partition function Z.k.i can readily be 
obtained by employing a trace invariance with respect to 
unitary transformation and using the diagonalized form of 
the elastic part of the bond Hamiltonian (IT) 



k,l 



^exp 

nt=0 



-pha)\nk + ^ 



2 sinh 



(9) 



Further, the magnetic part of the bond partition function 
Z,k.m can be replaced with a simpler equivalent expression 
provided by the generalized decoration-iteration transfor- 
mation ll23l|24i[25il26ll 



Zk,m = 2 exp 



= RoexpifiRio-l^o-l^), 



I2} 



cosh 



(10) 



which is of a general validity if and only if the mapping 
parameters Rq and Ri are given by the following formulas 



-(f) 



mi = - 



(11) 

(12) 



Substituting the elastic part of the bond partition func- 
tion (|9]l together with the generalized decoration-iteration 
transformation ( fTOl i into the formula (O, one readily gains 
an exact mapping relationship between the partition func- 
tion of the spin- 1/2 Ising model on a planar lattice with 
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the decorating atoms amenable to lattice vibrations and 
respectively, the partition function of the spin- 1/2 Ising 
model on a corresponding undecorated rigid lattice 



z 



exp(w)'^osh(f ) 



sinh2(^) 



(13) 



Thus, the partition function of the investigated model sys- 
tem can readily be obtained from the known exact results 
1 27 , 28 , 2^ 20, 31 ] for the partition function of the spin- 



1/2 Ising model on undecorated planar lattices with the 
effective interaction R\ given by Eq. (fT2] ). Besides, the 
critical behaviour of the spin- 1/2 Ising model on a deco- 
rated lattice with the vibrating decorating atoms can eas- 
ily be examined when comparing the effective nearest- 
neighbour coupling (fT2l) of the corresponding spin- 1/2 
Ising model on undecorated rigid lattice with its critical 
value (e.g. PARi] = 21n(l -i- V2) for a square lattice 
12711). It is also worth mentioning that the established 
mapping relationship (fTsT l between both partition func- 
tions permits a straightforward calculation of other ther- 
modynamic quantities such as the internal energy 



04 ^ nJL+HM ^ ^hojcoiM^^ 



(14) 



T 



16K 



(1 - 4eiM) + tanh (1 + 4eiM) 



J_ 

2 



\2 r 



Here, the former (latter) expression represents the lat- 
tice (magnetic) part of the overall internal energy and the 
quantity sim = (o'^i^ta^ labels the nearest-neighbour pair 
correlation function of the correspondin g spin- 1 / 2 Ising 
model on undecorated (rigid) lattice i29l bol bill . Note 
that the specific heat can be obtained from Eq. ST4\ by the 
use of well-known relation C = dUldT, but the final for- 
mula is too cumbersome to write it down here explicitly. 

Contrary to this, the spontaneous magnetization cannot 
be simply obtained from the mapping relation ( fTSl l be- 
tween the partition functions. The uniform and staggered 
magnetizations of the nodal Ising spins can be however 
derived by emp loying e xact mapping theorems developed 
by Barry et al. |32, 33, 34, 35t], which establish a precise 
mapping equivalence between the ensemble average (■ ■ ■) 
in the spin- 1/2 Ising model on a decorated lattice and the 
ensemble average (■ ■ ■)im in the corresponding spin-1/2 



Ising model on the undecorated rigid lattice for any func- 
tion depending just on the nodal Ising spins. Accordingly, 
the uniform and staggered magnetizations of the nodal 
Ising spins can be obtained from the mapping relations 

'"^ = ^('^H+'^ii) = ^('^h+'^L)im="^m(/5/?i >0K15) 

*A = \{^1-Q^k2) = \{^k-^k2)^= ^IMimi <0), (16) 

which connect both these quantities with the spontaneous 
(either uniform or staggered) magnetization of the corre- 
sponding^spin-1/2 Ising model on the undecorated rigid 
lattice IBafSTll . On the other hand, the spontaneous mag- 
netization of the decorating Ising spins niB can be related 
to the spontaneous magnetization of the nodal Ising spins 
mj with the aid of the generalized Callen-Suzuki identity 



niB 



tanh 



(^kl + ^k2^ 



- tanh 



(f)- 



Last, let us calculate the mean displacement and stan- 
dard deviation of the decorating atoms from their equilib- 
rium lattice positions. It is quite obvious from the coordi- 
nate transformation (|4]i that the mean displacement can be 
expressed in terms of two spin correlation functions and 
the mean displacement in a new coordinate system 



(17) 



It can be readily proved that both spin correlation func- 
tions between the decorating spin and its two nearest- 
neighbour nodal spins are identical 

(fil^li) = (fil^k2> = ^(1 +4eiM)tanh(yj (18) 

and the mean displacement in a shifted coordinate system 
equals zero 



2Mu 



(19) 



These results are taken to mean that the mean displace- 
ment also equals zero, i.e. (i^^) = 0, and the decorat- 
ing atoms oscillate symmetrically around their equilib- 
rium lattice positions. To bring an insight into a magni- 
tude of those lattice oscillations, let us also calculate the 
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standard deviation for the displacement of the decorating 
atoms from their equilibrium lattice positions using 



d = ^{xD - {h? = ^{(x'f) + (1 - 4eiM). (20) 

The mean of a square of the displacement of the decorat- 
ing atoms in a shifted coordinate system is given by 




and hence, the standard deviation for the displacement of 
the decorating atoms readily follows from the formula 

^=#'^°*(V)-^32^^^-^^™^- ^22) 

3. Results and discussion 

Now, let us proceed to an analysis of the most interest- 
ing results obtained in the foregoing section. Although all 
generic features of the investigated model system will be 
hereafter illustrated only on a particular example of the 
spin- 1/2 Ising model on a decorated square lattice with 
the coordination number q- A, the qualitatively same be- 
haviour may be expected for any other decorated loose- 
packed planar lattice as well. It is noteworthy, more- 
over, that the exchange and magnetoelastic constants can 
be considered without loss of generality positive, because 
the relevant change in a character of the exchange interac- 
tion does not fundamentally affect neither the critical be- 
haviour nor the behaviour of basic thermodynamic quan- 
tities as it merely causes a rather trivial change in a rel- 
ative orientation of the decorating and nodal Ising spins, 
respectively. Last but not least, let us introduce the follow- 
ing set of reduced parameters AjJ.KI J, hco/ J and k^TjJ, 
which measure a relative strength of the magnetoelastic 
constant, the spring stiffness constant, the phonon energy 
and the temperature, respectively. 

The typical ground-state phase diagram is depicted in 
Fig. |2]in the the JjK - AjJ plane. It turns out that the 
classical ferromagnetic phase (CFP) with a perfect par- 
allel alignment of all nodal as well as decorating spins 
constitutes the ground state when the magnetoelastic cou- 
pling is smaller than the boundary value A < A\, - y/8KJ. 




[K/J]-' 

Figure 2: The ground-state pha.se diagram in the J/K - A/ J plane. 

On the other hand, the more striking frustrated antiferro- 
magnetic phase (FAP) becomes the ground state when- 
ever the magnetoelastic coupling exceeds the boundary 
value A > Ab. The spin arrangement inherent to FAP 
can be characterized by a remarkable coexistence of a 
partial order and disorder, since the perfect antiferromag- 
netic (Neel) long-range order of the nodal spins is accom- 
panied with a complete disorder of the decorating spins. 
An origin of this rather unexpected and unusual phase 
apparently lies in the spin-phonon coupling induced spin 
frustration, which manifest itself in the competing charac- 
ter of the next-nearest-neighbour interaction between the 
nodal spins in the effective spin Hamiltonian ^ and leads 
to a frustration of the decorating spins. 

To bring an insight into how the phase transition be- 
tween CFP and FAP manifest itself in the lattice vi- 
brations of the decorating atoms. Fig. |3] displays the 
zero-temperature dependence of the standard deviation 
on a relative strength of the magnetoelastic constant. 
Even though the decorating atoms oscillate symmetrically 
around their equilibrium lattice positions both in CFP and 
FAP, it is quite evident from Fig. [3] that the magnitude 
of these lattice oscillations are much more robust in FAP 
than CFP. As a matter of fact, the standard deviation does 
not depend on the magnetoelastic constant within the re- 
gion corresponding to CFP where 

j,PP . (23) 

while it jumps towards a higher value at A = Ab and then 
increases monotonically with the magnetoelastic constant 
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within the region corresponding to FAP where 

From this point of view, our further analysis will be re- 
stricted only to a discussion of the particular cases with 
a sufficiently small value of the ratio A/K between the 
magnetoelastic and spring stiffness constants in order to 
keep a magnitude of lattice oscillations within the range 
of validity of the harmonic approximation. 

Now, let us proceed to a discussion of the critical be- 
haviour. In Fig. m the critical temperature is plotted 
against a relative strength of the magnetoelastic constant 
for several values of the spring stiffness constant. The 
critical temperature of CFP gradually decreases by in- 
creasing the magnetoelastic coupling until it completely 
vanishes at the ground-atate boundary between CFP and 
FAP. On the contrary, one may observe a rather steep in- 
crease of the critical temperature upon further increase of 
the magnetoelastic coupling in FAP, which can be mainly 
attributed to a quadratic dependence of the antiferromag- 
netic part of the effective interaction ( fTSl i on the magne- 
toelastic constant. It is worthwhile to recall, however, that 
the validity of our results for the critical temperatures of 
FAP is restricted only to a rather narrow region in a close 
vicinity of the ground-state phase boundary between CFP 
and FAP, where the standard deviation is small enough to 
ensure the validity of harmonic approximation. 

Next, let us turn our attention to temperature depen- 
dences of the spontaneous magnetization. Fig.|5ja) shows 



typical thermal variations of both spontaneous sublattice 
magnetizations in CFP for three different values of the 
magnetoelastic constant. The sublattice magnetization niA 
of the nodal Ising spins is always more resistent to thermal 
fluctuations than the sublattice magnetization mb of the 
decorating Ising spins even if both sublattice magnetiza- 
tions tend to zero with the same critical exponent /3p = 1/8 
from the standard Ising universality class. It is worthy to 
notice, moreover, that the difference in the thermal be- 
haviour of both sublattice magnetizations is the smaller, 
the stronger is a relative strength of the magnetoelastic 
constant. On the other hand, the relevant order parameter 
for partially ordered and partially disordered FAP repre- 
sents the staggered magnetization sa of the nodal Ising 
spins, while the decorating Ising spins remain disordered 
over the whole temperature range in FAP. The typical tem- 
perature variations of the staggered magnetization of the 
nodal Ising spins are depicted in Fig. Elb). It is worth 
mentioning that the thermal behaviour of the staggered 
magnetization in a close vicinity of the critical point is 
also governed by the critical exponent /3e = 1/8 from the 
standard Ising universality class. 

Last but not least, let us discuss the main features that 
are typical for temperature variations of the specific heat. 
Fig. |6] illustrates two typical thermal dependences of the 
specific heat for CFP, while Fig. [T] depicts another two 
typical temperature dependences of the specific heat for 
FAP. To gain a deeper insight into the overall behaviour, 
both lattice as well as magnetic part of the specific heat are 
plotted along with the total specific heat in those figures. 
It is quite obvious from Fig.|6ja) that the temperature vari- 
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Figure 5: (a) Thermal dependence of both spontaneous sublattice mag- 
netizations for the fixed value of the spring stiffness constant KjJ = 20 
and three different values of the magnetoelastic constant; (b) Thermal 
dependence of the spontaneous staggered magnetization of the nodal 
Ising spins for the fixed value of the spring stiffness constant KjJ = 20 
and two different values of the magnetoelastic constant. 



ations of the specific heat are initially governed almost ex- 
clusively by its magnetic contribution showing the /l-type 
singularity at the critical temperature, while the phonon 
contribution is just superimposed on the high-temperature 
tail of the magnetic heat capacity if considering relatively 
small values of the magnetoelastic constant. Another in- 
teresting non-monotonic thermal dependence of the spe- 
cific heat with the /l-type singularity superimposed on a 
round shoulder can be detected in Fig. |6|b), where the 
magnetic and phonon contributions to the overall heat ca- 
pacity are even better separated one from each other due 
to a relevant decrease in the critical temperature caused 
by an enhancement of the magnetoelastic coupling. This 
leads to an appearance of the round minimum in the high- 
temperature tail of the /l-type singularity before the spe- 
cific heat finally approaches its constant value C - INks 
in the limit of sufficiently high temperatures. 

Finally, Fig. |2la)-(b) illustrate typical thermal depen- 
dences of the specific heat in FAP, where a substantial 
shift in the critical temperature can be achieved by a rather 
small increase in a relative strength of the magnetoelas- 
tic constant. In this particular case, the lattice contri- 
bution to the overall heat capacity may prevail over the 
magnetic contribution even in the low-temperature region 
and hence, the /l-type singularity can be superimposed ei- 
ther on the ascending part (Fig.|7ta)) or the saturated part 
(Fig. |2|b)) of the lattice specific heat. Thus, the round 
shoulder in the low-temperature tail of the specific heat 
observed in Fig. [TJb) obviosly comes from the phonon 
rather than magnetic excitations. It is worthwhile to re- 
mark, moreover, that a few other types of dependences 
mostly with a rather well separated magnetic and lattice 
contributions can be achieved by varying the phonon en- 
ergy and the spring stiffness constant, however, our par- 
ticular attention was mainly concentrated on special cases 
with an obvious interplay between both contributions. 



4. Concluding remarks 



The present Letter deals with the magnetoelastic prop- 
erties of the spin- 1/2 Ising model with a spin-phonon cou- 
pling on decorated planar lattices partially prone to lattice 
vibrations, which are investigated by making use of the 
generalized decoration-iteration transformation and the 
harmonic approximation. Within the framework of this 
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Figure 6: Temperature dependences of the total specific heat, its mag- 
netic and lattice parts for the fixed value of the spring stiffness constant 
KIJ = 20, the phonon energy tliiilJ = 1 and two different values of the 
magnetoelastic constant: (a) AjJ = 5, (b) AjJ = 10. 
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approach, we have examined the ground-state and finite- 
temeprature phase diagrams, the mean and standard de- 
viation of the displacement of the decorating atoms from 
their equihbrium lattice positions, as well as, the temper- 
ature dependences of the spontaneous magnetization and 
the specific heat. 

It has been demonstrated that the spin-phonon inter- 
action gives rise to an efi'ective antiferromagnetic next- 
nearest-neighbour interaction between the nodal spins, 
which generally competes with the nearest-neighbour in- 
teraction between the decorating and nodal spins. As a 
result of this competition, one may encounter for suffi- 
ciently strong magnetoelastic couplings a quite intriguing 
partially ordered and partially disordered phase, where the 
antiferromagnetic long-range order of the nodal spins is 
accompanied with a complete disorder of the decorating 
spins. 

Before concluding, it is worthwhile to remark that 
the approach presented in this Letter can also be rather 
straightforwardly generalized to the mixed spin-(l/2, S) 
Ising model on decorated planar lattices with the deco- 
rating atoms of a quite general spin S > 1, where the 
spin-phonon coupling will produce an effective single-ion 
anisotropy and three-site four-spin interaction. Our future 
goal is to clarify the effect of the magnetoelastic coupling 
on a magnetic behaviour of these more general mixed- 
spin Ising models. 
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